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ABSTRACT 



Local and average heat transfer coefficients were deter- 
mined for a right circular cylinder in an oscillating flow. 
Spanwise platinum heater strips were used to heat the cylinder 
isothermally over the lower 180 degrees from the front to rear 
stagnation point. The four inch diameter cylinder was posi- 
tioned both normal to and at 45 degrees to the flow direction. 

Data was gathered for diameter Reynolds numbers from 
100,000 to 300,000. Large amplitude oscillations were imposed 
upon the mean flow using a rotating shutter arrangement. Fre- 
quencies of 0, 5, 10, 22, 50, 100 and 126 Hertz were investigated. 

For normal flows, local heat transfer coefficients in the 
wake and average Nusselt numbers were enhanced above values 
for steady flow, for runs having large amplitude. In the 45° 
flows, no significant change in heat transfer was noted with 
oscillating flow. 
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I. 



INTRODUCTION 



A. NORMAL CYLINDER 

Forced convective heat transfer from a cylinder in uniform 
flow has been a topic of interest to researchers since the turn 
of the century. In 1914 Professor L.V. King presented his 
relationship to describe the heat transfer from a cylinder in 
laminar crossflow as a function of air velocity: 

Nu = A.. + B-.Nt, 0,5 

1 1 Re 

Although slightly modified by Collis and Williams in 1959, 
it still serves today and is the basic equation of hot wire 
anemometery. Numerous other correlations are available and 
are a standard part of the heat transfer literature [Ref. 1]. 

A topic of increasingly greater interest to researchers 
and designers is the effect of an oscillating flow on forced 
convection. Gas turbine designers are intensely interested 
in upgrading the performance of gas turbine blades through 
materials research and blade cooling. Gas turbine blades 
experience a hot, oscillating flow environment as they pass 
through the wakes of upstream blades. 

Heat exchange designers involved in waste heat recovery 
from diesel engine exhausts are also concerned with the effect 
of flow oscillations on the exchange of heat to and from 
cylinders . 
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Little data has been made available to date for a cylinder 
in oscillating flow. Data is available for heat transfer from 
a flat plate in oscillating flow [Refs. 1, 2 and 3] and for 

an airfoil in oscillating flow [Ref. 4]. Data for a cylinder 
in other than uniform flow is generally for a vibrating cylinder 
with or without a net flow [Ref. 5] or for a cylinder with 
screen induced turbulence of less than about ten percent [Refs. 

6, 7, 8 and 9]. This study concerns itself with the effect of 

oscillations on the order of fifty percent of the freestream 
velocity on local heat transfer rates. 

Yawed cylinders are of interest to aerodynamicis ts because 
of the correlation of a yawed cylinder to a swept back aircraft 
wing. Researchers and designers may also encounter yawed cylinders 
in heat exchangers and other applications. 

No information is currently available on the effect of 
flow oscillations on a yawed cylinder. Viteri [Ref. 10] deter- 
mined the average heat transfer coefficients of yawed cylinders 
in uniform flow, and Kraabel [Ref. 11] determined the local 
heat transfer coefficients around yawed cylinders in uniform 
flow. These papers provide a basis for validation of the yawed 
cylinder results reported here. 

It has been demonstrated by Miller [Ref. 2] that augmented 
convective heat transfer can be attributed to the reduction of 
laminar to turbulent transition Reynolds numbers. Despard 
[Ref. 12] showed that freestream oscillations promote laminar 
boundary layer separation in the presence of an adverse pressure 
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gradient. Banning [Ref. 13] produced data that revealed that 
the starting vortex associated with periodic separation from 
an airfoil in oscillating flow produced an unexpectedly strong 
boundary layer reattachment which markedly improved the aero- 
dynamic performance of the foil. 

In summary, the effects of early laminar to turbulent boun- 
dary layer transition and the destabilized boundary layer pro- 
duced slightly augmented heat transfer while periodic vortices 
shed from the separation point can cause periodic reattachment 
of the hydrodynamic boundary layer. This can be expected to 
greatly enhance the heat transfer in the afterbody region. 

The purpose of this investigation is to determine the effect 
of large flow oscillations on the local heat transfer coeffi- 
cients around both normal and yawed cylinders. Additionally it 
is desired to develop a correlation formula for the average 
Nusselt number in oscillating flow that includes the effect 
of yaw angle, Reynolds number, and flow frequency and amplitude. 
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II. DESCRIPTION OF APPARATUS 



A. MODEL 

The model used in the investigation was a four inch dia- 
meter cylinder. It was fabricated from one half inch wall 
thickness fiberglass-phenolic tubing two feet long. The 12 
inch mid-span section was instrumented with 16 platinum heater/ 
resistance thermometers oriented parallel to the axis of the 
cylinder and spaced 12 degrees apart from center-to-cen t er . A 
sketch of the model shox^ing the orientation of the heater strips 
is shown in Figure 1. 

The platinum strips used were 12 inches long by one-quarter 
of an inch wide by one- thousand th of an inch thick. Each strip 
subtended an arc of 7.16 degrees of the circumference of the 
cylinder. The strips were spot welded to brass studs that were 
machined flush with the cylinder on the outside and protruded 
into the interior of the model for attachment of the voltage 
sensing and current carrying wires. Figure 2 shows a cross 
section of the cylinder wall through the studs and strips. 

Voltage drop across the platinum strip was measured using the 
voltage sensing wires attached to each stud. Platinum was 
used to both heat the cylinder and determine the surface tem- 
perature of the cylinder because of its large linear resistivity- 
temperature relationship. The resistance of the strips was 
determined from the voltage drop across the strip which was 
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then to be used to find the temperature. The amount of heat 
convected at each location was proportional to the power dis- 
sipated by the strip which was equal to the product of the cur- 
rent squared and the strip resistance. 

The entire outside of the cylinder, including the platinum 
strips, was covered with a thin layer of hard epoxy. The hard 
epoxy layer was machined and polished on a lathe to reduce the 
thickness over the strips to less than one-hundredth of an inch 
and to produce a smooth heat transfer surface. A photograph 
of the completed model showing the mounting stud and bracket 
and electrical leads is shown in Figure 3. 

Two copper-cons tan tan thermocouples were mounted out. the 
inside surface of the cylinder at 90 degrees and 270 degrees 
from the front stagnation point. These thermocouples were used 
to monitor the interior wall temperature of the cylinder in order 
to estimate the heat loss to the inside of the cylinder. 

B. WIND TUNNEL 

The experiment was carried out in the oscillating flow wind 
tunnel of the Naval Postgraduate School. Photographs of the 
facility are shown in Figures 4 and 5. The oscillating flow 
wind tunnel is a low speed, open circuit wind tunnel with the 
inlet bell extending out of the building. The tunnel inlet 
bell is eight feet square and the test section is two feet 
square, providing a 16.1 contraction ratio. Three high solidity 
screens located in the inlet section upstream of the nozzle 
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produce freestream turbulence intensities of less than 0.5 
percent of the freestream velocity for uniform flow. 

The wind tunnel drive consists of two axial flow fans in 
series, each of which has an internal, 100 horsepower, direct 
connected, 1750 RPM motor. The fan blades are internally adjust- 
able through a pitch range of 25 to 55 degrees, providing a 
wide operating base. A set of variable inlet vanes, located 
immediately upstream of each fan, are externally adjusted to 
provide control of test section velocity. These vanes, of radial 
configuration, preswirl the air in the direction of fan rotation 
to reduce fan capacity. The total range of freestream velo- 
cities possible is from ten to 250 feet per second. 

A plan view sketch of the wind tunnel showing the location 
of the fans and inlet vanes in relation to the test section is 
shown in Figure 6. 

Two fundamental methods are available to superimpose oscil- 
lating flow, on a mean velocity. Both Nickerson [14] and 
Hori [15] introduced oscillations by oscillating their models 
in a steady flow environment. This method severely restricts 
the range of attainable frequencies because of mechanical com- 
plications and instrumentation difficulties. 

The other approach is to provide an oscillating flow environ- 
ment over a stationary model. Hill [16] used a sliding shutter 
to impose oscillations on the freestream but was limited to 
low frequencies because of mechanical difficulties. Feiler 
and Yeager [17] used a siren mounted upstream of an eight inch 
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diameter wind tunnel to expose small models to frequencies 
ranging from 34 to 300 Hz. with root-mean-squared amplitudes 
of up to 65 percent of the mean velocity. The most successful 
method of obtaining oscillating flows with a large range of 
frequencies, amplitudes and flow velocities in a moderate sized 
tunnel is that employed by Karlson [18] and later by Miller 
[19]. A rotating shutter valve, immediately downstream of the 
test section is used to superimpose a periodic variation of 
velocity on the mean flow. The method used in this investiga- 
tion is identical to that used by Miller. 

The shutters consist of four horizontal, equally spaced, 
shafts mounted in a plane just downstream of the test section. 
The shafts are slotted to receive blades of various widths, 
forming a set of four butterfly valves spanning the wind tunnel. 
Figure 7 is a picture of the shutter valve assembly. The 
shutters are driven synchronously by a pulley and belt arrange- 
ment and a vari-drive motor and gearcase, and can be varied 
from 0.1 - 250 Hz. The amplitude of the oscillation is a func- 
tion of the width of shutter blades, frequency of rotation of 
the shutters and the freestream velocity. Five sets of blades 
are available producing tunnel blockages of 33, 50, 66, 83 

and 100 percent of the flow area in the fully closed position. 

In this investigation blades providing 66 and 100 percent 
blockage were used resulting in amplitudes of 5 to 69 percent 
of the mean freestream velocity. 
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The test section of the tunnel is fabricated from two-two 



inch thick, 25 inch wide by 18 feet long pieces of aluminum 
as the upper and lower walls and three-two inch thick stress 
relieved lucite panels on each side as side walls. The front 
panels are hinged at the top and can be hydraulically opened 
and closed for access to the test section. 

For this investigation the rear center, lucite panel has 
been replaced by a two inch thick plywood panel to facilitate 
mounting the model. The model was mounted by a chordwise 
extension held in a rotatable clamp mounted flush in the center 
of the plywood panel. The other end of the model is supported 
by a 3/8 inch thick "Y" shaped, faired bracket, that is affixed 
to the tunnel ceiling and floor, and rests against the front 
lucite panel. The model is shown in place with the M Y" bracket 
clearly visible in Figure 8. Electrical leads were taken out 
through the plywood back panel where they were connected to the 
power supplies and the voltage drop sensing circuits. The 
model was yawed by affixing extension pieces to each end that 
were angled and sized to completely span the tunnel at each 
angle of yaw. 

Power to heat the strips was controlled by sixteen variacs 
that powered 16 individual DC power supplies (see working cir- 
cuit, Figure 9). A precision shunt (R=0.001 ) was placed in 
series with each strip and the voltage drop across the known 
resistance was used to determine the current. The voltage 
drop across each strip was measured using the potential sensing 
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leads that were attached to the studs on the inside of the model. 



Instantaneous velocity measurements for the purpose of free- 
stream turbulence determination in uniform flows and amplitude 
and frequency measurement in oscillating flows were made using 
a linearized constant temperature hot-wire anemometer described 
in Reference 20. An oscilloscope and a Ballantine True R.M.S. 
meter were used to monitor hot-wire signals. Mean velocity of 
the free stream was measured using a Prandtl Probe and a water- 
filled micro manometer. 

Figure 10 shows a pair a typical velocity wave forms. The 
vertical axis is velocity and the horizontal axis is time. The 
roughness of the upper tracing is typical of low amplitude oscil- 
lations which are significantly more difficult to produce in 
true sinusodial form. 

The freestream temperature and the model internal wall 
temperature were measured using copper-cons tan tan thermocouples 
and a Newport digital temperature indicator. 

The frequency of oscillation was also monitored with an 
optical chopper driven by a shutter valve shaft and a digital 
counter . 
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III. EXPERIMENTAL PROCEDURE 



A. CALIBRATION OF MODEL 

The platinum strips were calibrated as resistance thermo- 
meters by determining their electrical resistance as a function 
of temperature in an electric oven. Seven separate tempera- 
tures were used, evenly spaced between ambient (60°F) and 
200°F. Although the platinum h ea t er- the rmome t e r strips were 
extremely pure and are expected to have a linear resistance- 
temperature characteristic, calibration was carried out to 
determine the exact slope and intercept of each strip. 

The calibration circuit employed (see Figure 9) was powered 
by a low impedance DC voltmeter calibration standard. This 
proved to be an extremely precise, drift free power supply. The 
voltage drop across the precision resistance (1.0 + 0.15%) was 

recorded as was the voltage drop across the heater strip itself. 
The current was then computed from the voltage drop across the 
precision resistance divided by one ohm, and the resistance of 
the strip for that temperature was computed from the voltage 
drop across the strip divided by the current. Temperatures 
in the oven were monitored with a Newport Laboratories Digital 
Temperature Indicator and a cop p e r- cons t an t an thermocouple. 

After determining the resistances as a function of the 
temperature a straight line was fitted to the data using linear 
regression. The resulting equations were used throughout the 
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subsequent experimentation and data reduction to determine sur- 
face temperatures. 

B. ADJUSTMENT OF SURFACE HEATERS 

Before each series of runs the heater strips were energized 
and the wind tunnel was run with a low velocity uniform flow 
in order to allow the model substrate to come to thermal equili- 
brium. Once this was accomplished the shutters were turned on 
and adjusted to the desired frequency and the velocity was 
adjusted to the magnitude corresponding to the Reynolds number 
desired. At this point the model was again allowed to come to 
thermal equilibrium. 

Once this initial thermal equilibrium was achieved the 
voltage drop and current through each strip was recorded and 
inserted into a predictor-corrector program (see Appendix F), 
This program calculated a new required voltage drop based on 
desired strip temperature, current strip temperature, and the 
known r esi s tance- t empe rature calibration curve for that strip; 
the strip was then adjusted to the new desired voltage drop 
through adjustment of its power supply. Following adjustment 
of the strips the model was again allowed to return to thermal 
equilibrium. 

These adjustments were performed as many times as necessary 
until all the strips were at the same temperature, ensuring an 
isothermal heat transfer surface that minimized heat conduction 
between strips . 
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The temperature for a particular run was determined from 
a sampling of the strip temperatures before the first set of 
adjustments was made and was selected to minimize the number 
of adjustments required. In every case however the temperature 
was chosen to be over 110°F in order to maintain at least a 
50° temperature difference between the strip and the ambient 
tempe ra ture . 

C. FREESTREAM CONDITIONS 

A series of runs was made with both the 45° yawed model 
and the normal model covering a wide range of operating condi- 
tions. Each model was tested with both the six inch blades 
(100 percent tunnel blockage) and the four inch blades (66 per- 
cent tunnel blockage), resulting in oscillation amplitudes from 
five to 69 percent of freestream mean velocity. Each combina- 
tion of frequency and amplitude was run at two freestream velo- 
cities corresponding to nominal diametral Reynolds numbers of 
300,000 and 150,000. Frequencies investigated were 0, 5, 10, 

22, 50, 100 and 125 Hz. 

D. DATA REDUCTION 

The data reduction technique is described in Appendix A 
and a sample calculation is presented in Appendix B. A detailed 
energy balance is necessary in order to determine what portion 
of the total Joulian heat generated by platinum resistance 
heater is transferred from the surface by convection. Part of 
the energy was conducted through the thin epoxy coating and 
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convected directly to the f reestream or was radiated into the 
environment. The balance of the energy was conducted into the 
substrate of the cylinder where some was conducted and con- 
vected into the center of the cylinder, some was lost to the 
ends of the cylinder, and some was conducted to the gaps between 
the strips where it was convected to the freestream or was radia- 
ted to the environment. 

The energy lost from the substrate side of the surface 
heaters is hard to determine analytically since it consists of 
both conduction to the intraheater gaps and subsequent loss by 
convection and radiation as well as conduction to the inner 
surface of the model. Therefore a two-dimensional Teledeltos 
paper model was constructed to determine the relationship 
between the thermal resistances to the inner and outer surfaces. 
This procedure and the results are discussed in Appendix D. 

The results of the Teledeltos paper experiment were placed in 
equation form and solved interatively for the external covec- 
tion heat transfer coefficient, h, for each strip. 

The other major loss of energy from the system is by radia- 
tion to the environment. These losses can be appreciable and 
an accounting for this loss was developed as part of the data 
reduction program. Details are described in Appendix C. 
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IV. RESULTS AND DISCUSSION 



A. STEADY FLOW 

Six steady flow runs were made, two with the model normal 

to the flow at N_ = 150,000 and 300,000 and four runs with 

Re 

the model at 45° to the flow at N^ = 35,000, 100,000, 150,000, 

and 300,000. The results of these runs are shown graphically 
in Figures 13 through 18. Results of the model at 90° (normal 
to the flow) at N R ^ = 150,000 are compared in Figure 19 with 

the theory of Froessling [Ref. 24] and the experimental 
results of Schmidt [Ref. 25] for N R ^ =• 170,000. The discrepancy 

between theory and both of the experimental results is expected 
and is a consequence of the effect of freestream turbulence 
near the stagnation point [Ref. 26]. 

Comparison of the results of the present experiment with 
those of Schmidt indicates two interesting points. The two 
agree within 20 percent up to an angle of 150° from the 
stagnation point, beyond which they diverge markedly, differing 
by 70 percent at 180°. This sharp increase in heat transfer 
above that noted by Schmidt is typical of all the results 
observed in the present investigation and shows no correlation 
with frequency number or Reynolds number. Agreement within 
20 percent is good however and this lends a high degree of 
confidence in the experimental techniques and results. 
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No local heat transfer data is available for the steady 
flow case with 45° yaw. However, the data of Viteri [Ref. 10] 
can be interpolated to arrive at an average Nusselt number of 
128.3 at a Reynolds number of 31,900. At a Reynolds number of 
35,630 the present experiment produced an average Nusselt number 
of 122.1. This close agreement lends confidence in at least 
the average Nusselt Numbers obtained in the 45° yawed case. 

B. NONSTEADY FLOW 

Two thirds of the nonsteady freestream conditions did not 

produce substantial alteration in heat transfer when compared 

to the steady flow case. Figures 20 through 25 depict the 

cases with the 90° model that show little or no altered heat 

o 

transfer, and Figures 26 through 39 depict those of the 45 
model. In each of these figures the steady flow data has been 
plotted along with the nonsteady data for comparison. 

Twelve of the 38 freestream conditions investigated pro- 
duced substantially altered local heat transfer rates. These 
results are shown graphically in Figures 40 through 49 for the 
90° model, and Figures 50 and 51 for the 45° model. The 
changes in local Nusselt number took place almost exclusively 
in the afterbody region of the model, and resulted from altered 
separation and flow patterns in the wake of the cylinder. Fig- 
ure 40 is a typical example of these results and demonstrates 
enhanced heat transfer in the wake area of the 90° model. The 
data reported in Figures 48 through 51 exhibit diminished heat 
transfer rates. Figure 50 is typical of this class of data. 
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In order to simplify the analysis of the results, the 
average Nusselt number was calculated for each freestream con- 
dition (see Table I) . Figures 52 and 53 report average Nusselt 
numbers versus Reynolds number for the 90° model and the 45° 
model respectively. In the 90° case (Figure 52) a curve has 
been drawn representing Nu = N^ and the Amplitude and 

Storuhal numbers have been annotated for each data point. Of 
interest is the fact that data with higher Amplitude numbers 
are those which tend to exhibit augmented average heat transfer 
rates. Moreover, higher Amplitude numbers seem to correspond 
to greater heat transfer augmentation. Although the data is 
sparse, a pattern of larger Nusselt numbers for larger Ampli- 
tude numbers at the same Reynolds number seems to exist. Three 
apparent anomalies exist associated with the data from run num- 
bers 15, 17 and 20. It is possible that an error in recording 

the amplitude data was made and that the Amplitude numbers for 
runs 17 and 20 were reversed. If this were true then both 
runs 17 and 20 would conform to the apparent pattern. When the 
Frequency numbers, Nusselt numbers and Reynolds numbers of runs 
15 and 17 are compared, it is seen that they are very similar 
but that the Amplitude numbers are considerably different. If 
the previously noted pattern for higher Nusselt number for 
higher Amplitude number is true then the Amplitude number of 
run 15 is suspect and may be too low. 
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Average Nusselt number data taken for the 45° model is 
shown in Figure 53. No apparent pattern of heat transfer altera- 
tion is evident in this plot even though Amplitude numbers as 
high as 0.44 and 0.65 appear in the data. 

C. NUSSELT NUMBER CORRELATION 

The data obtained in this investigation is too sparse to 
allow a correlation between heat transfer rate and freestream 
conditions to emerge. Only a possible pattern is suggested by 
the results and more data covering a wide range of Reynolds 
numbers and Amplitude numbers is necessary to produce such a 
correlation . 
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V. CONCLUSIONS 



Although there is insufficient data to draw any definite 
conclusions or to derive a correlation, several trends may be 
observed. ,In the 90° case, significant increases in the local 
Nusselt numbers are seen in the wake region and are due to 
changes in the separation and wake structure. Moreover, aug- 
mentation of local Nusselt numbers seems to increase with 
increasing frequency as suggested by the data shown in Figures 
40 through 47. Increases in average Nusselt numbers is appar- 
ently related to the Amplitude numbers as discussed above and 
shown in Figure 52. With the exception of the three points 
associated with runs 15, 17 and 20, discussed in Chapter IV, 

the pattern exhibited is for higher average Nusselt numbers to 
be associated with greater amplitudes. 

In the case of the 45° yawed cylinder, no general patterns 
are discernable and it appears that local Nusselt numbers 
exhibit no significant change over a wide range of frequencies 
and amplitudes. 

This investigation has served to initiate a study of heat 
transfer from a cylinder in oscillating flow. Additionally 
it has validated the techniques to be used for further study. 
More extensive data collection will be required before defini- 
tive conclusions can be drawn or correlations made. 
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APPENDIX A 



DATA REDUCTION 

The data was reduced to the form shown in Table 2 by a 
specially written Fortran computer program on an IBM 360 com- 
puter at the Naval Postgraduate School. 

The raw data consisted of pitot-static tube £P (inches 

HO), freestream temperature, T , (°F) , internal model tempera- 

Z r 

tures (heater side and "other" side, °F), frequency of oscil- 
lation (Hz), amplitude ratio (%) and the voltage drop across 
(volts) and current passing through (amperes) each strip. The 
total power dissipated by each strip is simply: 

P = 3.4149 x E x I 

where 3.4149 is a conversion factor from watts to BTU/Hr. 

From the measured current (I) and voltage drop (E) each 
individual strip temperature is determined by comparing its 
resistance to its known temperature versus resistance charac- 
teristics. The resistance of a strip is: 




From the linear temperature-resistance relationship obtained 
during calibration, the strip temperature can be obtained: 
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where C and S are the known ordinate intercept and slope 



obtained from the individual resistance versus temperature 
calibration curve. 

The other unknown temperature required by the data reduc- 
tion program is the inner wall temperature, The inner wall 

temperatures are known at two points. Thermocouples are mounted 
on the wall at locations 90° and 270° from the stagnation point 
at mid-span (see Figure 2). The inner wall temperature is 
assumed to vary linearly between these two measured values. 

With this information h, the heat transfer coefficient may 
be calculated. The following formula (developed in Appendix 
D) is used : 



■ 

T - T_ 
s F 


Rt, 

1 + 


T - T 
S I 


1 

M 

1 


R^ 






R_ 


D 


0 

m - 


D 


D 

«■ 



h is contained in the expressions for R^ , R q /R d , and 
R^/Rp as is shown in Appendix D. This equation was solved 
iteratively to obtain h. 

With the local h for each strip, the local Nusselt number 
is calculated from: 



Nu 



hD 

k 
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APPENDIX B 



SAMPLE CALCULATION 



For the purpose of this sample calculation strip Nr 6 from 
data run Nr 20 was selected. 



The information recorded is: 
Voltage drop 
Current 

Pitot-Static Pressure 

Amb i en t Temp er a t ur e 

Inside Temperature : 
Heater Side 
Other Side 

Frequency 

Amplitude Ratio 

Blades 

Model 

Barometric Pressure 
Freestream Mean Velocity 



V 



■[ 



2g«W (£ff£fE) - 1) 



] 



1/2 



E = 1.62 Volts 
I = 7.47 Amps 



HW = 4.16 inches H 9 0 



TA = 64.0 F 



TIHTR = 105.0 F 
TIOTH = 7 7 . 0°F 



Freq = 50.0 CPS 

N. = 0.31 
A 

Medium (4") 

90° to stream 
30.04 inches Hg 



62.4 



2(32.2 •£-!•) (4.16 in H O) 

S 2 2 12 - in | 0.07654 



lb f 
ft 3 



lb f 
ft 3 



- 1 



1/2 



= 134.4 



ft 

sec 
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Freestream Reynolds number: 



j = ™ 
Re v 



134.8 — (0.333 ft) 
s 

0.000166 f 1 2 / s 



= 269,800 



Temperature of the strip 



T = 



!-«= 



E = 1.622 Volts 
I 7.47 Amp s 



= 0.2171 ohms 



C = 0.173955807 



S = 0.000390625 



Therefore , 



T = 



0.2171-0, 173955807 
0.000390625 



= 110.5 F 



Total Heat Dissipated by the strip: 



q s = El = (3.4149 



Btu 



Vo 1 t-Amp-hr 



>(1.622 Volts) (7.47 Amps) 



= 41.38 Btu 



Also : 



q s R 



D L 



< T S - V' 1 + R^> + < T S 



V<iq> 



r d - 



k + ih 

] 

kh T A 



r t 



a l + a 2 h T 
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R d b l h T + b 2 h T + b 3 h T 



After substitution, the equation may be solved for . 
(the constants a^ , a^, b^, b ? and b^ are given in Appendix D) 

By iteration: 



h = 36.9 



B tu 



hr - ft 2 - ° F 



The radiant heat transfer coefficient 



h r - e<r(T s + t a )(t s 2 + x A 2 ) 



Assuming is the same as the fluid temperature: 

h r = (0. 73) (0.1718 x 10~ 8 ) (570.5 + 524.0) 

x (570. 5 2 + 524. 0 2 ) 



= 0.824 



B tu 



hr - ft 2 - °F 



Therefore, the local convective heat transfer coefficient is 



h = b - h 
T r 



= 36.9 - 0.824 
Btu 



= 36.1 



hr - ft 2 - °F 
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The local Nusselt number: 



hd _ (36.1) (0.333) 
k ~ 0.014738 



= 815.9 
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APPENDIX C 



RADIATION 

For a uniformly convex gray body with emis s ivi ty , at a 
temperature T, radiating into an enclosure of temperature T — , 
the net radiant heat exchange is: 

q r = <Te A(T 4 - T^ 4 ) 

where ^ is the Stef an-Boltzmann constant, and A is the surface 
area of the emitter. 

The above equation can be put in the same form as the con- 
vection equation: 



q = h A (T - T ) 

r r ® 

where h is the effective radiation heat transfer coefficient 
r 

given by: 

h r = (T + T te ) (T 2 + T^ 2 ) 



Taking £ = 0.73 for epoxy coating over aluminum [Ref. 22] 
and the following typical temperatures: 

T = 120°F = 580°R 
T # - 60°F = 5 20°R 
Gives : 



h r = 0.73 (0.1714) (10 8 )(580 + 520)(580 2 + 520 2 ) 



= 0.835 



Btu 

hr • f t 2 • F 
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This value falls within two to 10 percent of the total 
heat transfer coefficient and therefore cannot be neglected. 
Therefore, the convective heat transfer coefficient, h, is 
obtained by subtracting the radiant heat transfer coefficient, 
h , from the experimentally derived total heat transfer coef- 
ficient , that is , 



h = 
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APPENDIX D 



CONDUCTION LOSSES 

Although the thermal resistance of the fiberglass substrate 
is very high compared to the thermal resistance between the 
strip and the fluid, some heat is lost to the interior of the 
model. This heat loss must be deducted from the total heat 
input in order to arrive at the actual heat transferred to the 
air. Heat loss due to conduction in the longitudinal axis direc- 
tion and radially through the small brass studs was neglected. 

A Teledeltos paper model was used to estimate the heat loss 
to the interior of the model. Symmetry was used to reduce the 
size of the required model (see Figure 11). 

Teledeltos paper is a graphite covered paper having a nomi- 
nal resistivity of 1000 ohms per inch square and is highly 
isotropic. Teledeltos paper is extremely useful for two- 
dimensional heat conduction modeling in cases such as this 
which have boundary conditions and a physical shape that are 
inconvenient for an exact solution. 

The heater strip itself is simulated with a strip of silver 
conductive paint that is placed on the paper in exactly the 
same position as the platinum strip is on the actual model. 

The inside of the model segment is considered isothermal so 
the analogous surface of the Teledeltos paper model was 
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held at uniform potential with a conductive strip of silver 
paint. The left and right hand adiabatic boundaries, between 
strips, are easily modeled by leaving them bare. The convective 
boundary at the gap between strips is modeled by allowing the 
Teledeltos paper to extend beyond the actual cylinder surface 
for a distance to simulate the convection resistance. The 
required distance is determined by: 




where k is the conductivity of the actual substrate, h, is the 
convective coefficient being modeled and S is the scale factor, 
in this case, S = 48. Additionally the excess length was slit 
into eight equal width strips lengthwise to force the flux 
lines parallel to each other and perpendicular to the surface. 
In order to model h T s of 10, 20, 50 and 100 Btu/hr • f t^ • F , L f s 

of 14.4, 7.2, 2.88 and 1.44 inches respectively were used. 

A strip temperature of 120°F with an inner temperature 
of 100°F and a fluid temperature of 60°F were used. These 
temperatures were scaled to appropriate voltages that would 
avoid excessive Joulian heating of the Teledeltos paper. 

In operation the voltages were measured at 17 equally 
spaced intervals along the line corresponding to the cylinder 
surface and at five equally spaced points along a line cor- 
responding to 0.08 inches from the inner surface. This allows 
the determination of the ratio of heat loss to the inside of 
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the model, compared to the heat loss indirectly through the 
air (via the substrate). 

Considering the lumped electrical resistance network 
analogous to the heat transfer problem (see Figure 12) it is 
noted that the power input to the strip is equal to the sum of 
the losses to the air, to the inside of the cylinder to the 
gap between the strips, or: 

“s ■ q D + «0 + ’I 



since in general: 



4 T 




this becomes : 




Where T , Tp and are the temperatures of the strip, the 
fluid and the inside respectively and , R Q and R-j. are the 
thermal resistances from the strip to the air, from the strip 
through the gap to the air (R Q = R a + R b ) and through the 
substrate to the inside of the cylinder respectively (see 
Figure 12 ) . 

This equation can be manipulated into the following form: 
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The resistance to heat loss from the strip directly to the 
air, , is due to the thin layer of epoxy on the strip and 
the convection film resistance. The conduction resistance 
due to the epoxy is equal to -£/kA where is the thickness of 
the epoxy layer, k is the thermal conductivity of the epoxy 
and A is the surface area of the strip. The convective resis- 
tance is equal to 1/hA, where h is the convection heat transfer 
coefficient. Summing these we obtain: 

R = k ±Jh 

D khA 

The quantities R^/R^ and R^/R^ are derived experimentally 
from the Teledeltos paper analogue and are in the form of 
polynomial equations in h: 

R o 

^ = 1.546744 + 0.0320672 h 

R D 

R I 2 3 

^ = 0 . 106105h + 0 . 0047559ft + 0.0000079h 

R D 

These resistances were used in the heat transfer equation 
to solve for h (see Appendix B) . 
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APPENDIX E 



UNCERTAINTY ANALYSIS 



For the purpose of this analysis, the method of Kline and 
McClintock [Ref. 23] will be used, assuming that all variables 
have a normal Gaussian distribution. Possible dimensional error 
in strip size, placement and other dimensions were considered 
so small as to be negligible. The tabulated values of k and 
e were assumed to be correct. All temperatures measured with 
thermocouples were assumed to have an uncertainty, due to 
readout error, of 0.1 degree F. 

The strip temperature uncertainties derive from the 
+0.003 mV uncertainty in reading the digital voltmeter used in 
calibration and the +0.002 V and +0.01 amp uncertainty in 
reading the meters during data runs. 

The uncertainty in measured resistances during calibration 
is 



= +0.00031 ohms 

The resistances obtained in calibration are used to obtain 
a straight line equation for temperature vs resistance. 

R = S-T + C 




2 



1/2 
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The uncertainty of C is assumed to be of the same order as 
R. The uncertainty of S is: 



A S = S 



( 4 R )2 + 



-1 1/2 



= +0.00000064 



Since the strip temperature during a data run is given by 



- C 



T = 



The uncertainty of T is given by 



A T = T 
S S 



( ^|) 2 + (^ I ) 2 + ( r | S ) 2 



- 1/2 



+ T 



(^) 2 + (^) 2 



1 1/2 



= +0.54 F 



The values of R^/R^ were obtained from a Teledeltos paper 
model. Uncertainties in measuring a potential difference between 
any two points on the Teledeltos paper was +0.02 volts. This 
corresponds to a uncertainty in temperature of +0.2°. For 
Rq / Rq is given by: 

R 0 _ APT 

R d sadt 
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where ADT is the direct convection surface area times (T - T_) 

o r 

and SADT is the indirect convection area times (T local - T„) . 

r 

(T local is determined from the Teledeltos paper model.) 

The uncertainty of SADT is given by: 

SADT = SADT ( ^p -~) 

= + 0.0011 



Similarly for ADT 



A ADT = ADT (^|) 



= + 0.0021 



Now 



R T 



( AADT ) 2 + (ASADTp 



n 1/2 



ADT' 



SADT' 



= + 0.0088 



Rj/Rp is given by: 

R I _ ADT h ^S ~ T I ) 

R d IADT k - T p ) 

The uncertainty of IADT is given by: 

A IADT = IADT ) 

= + 0.032 
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Thus the uncertainty of R^/R^ is given by: 



A 





-4 ADT x 2 ,AlADT 2 

ADT } + (_ IADT ) + 




+ 



\ rp rp / 

S F 



1/2 



= +0.024 



q is given by : 



q = El 



The uncertainty in q is given by: 



Aq = q 



1/2 



(^ f ) 2 ♦ c ^) 2 



= +0.052 



Btu 

hr 



Because of the complexity of the formula used to determine h, 
h was found by iteration of the following implicit equation: 



q 



hkA 

k+hjL 




T F^ (1 + A1 + A2h 



) + 



(T 
V S 



T x ) ( 



Blh + B2h 2 + B 3h 3 



The uncertainty of h was determined by: 



AT 



Ah = 



S2h 



2T. 



2 AT 2 4T I2h 2 

-) + ( F S~) + ( H-) 



2T, 



2T. 



n 1/2 



+ (A(^) —111 ) + 2h 

( ( R n } 2 (R q /R d ) ) C ( R n } 2 ( R T / R^ ) 



V D' 
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The partial derivatives were evaluated numerically and were 
found to be : 




= -0.9137 




= 0.0698 



2h 

ZCRj/Rjj) 



0.01 




= 1.375 



2h 

2 (Rq /R D ) 



4 . 7619 



therefore : 



Ah = +0.51 



B tu 

Hr . f t 2 • F 



The mean velocity is given by: 



V 



M 




water 

air 




1/2 



The only parameter here that is open to uncertainty (due to 
the assumption that property values are correct) is the H, 
or height of water column in the micromanometer. Since 
AH = +0.02 inches: 



„ .. M ,AH N 

AV M 2 ^ H 



= +0.324 ft/s 



46 



since : 



N 



Re 



VD 

V 



AL = N„ 
Re Re 



V 

( A_M) 

V V ’ 



M 



= + 635 . 



Nusselt number is by definition: 



Nu 



hD 

k 



therefore : 



A Nu 



Nu 




= + 11.3 
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APPENDIX F 



PREDICTOR - CORRECTOR EQUATION 

In order to simplify and shorten the process of adjusting 
strip voltages to obtain an isothermal surface as described 
in Chapter 3, a predictor-corrector program was developed and 
implemented on a TI-59 with PC-100C thermal printer. The pro- 
gram listing is included at the end of this Appendix. 

The program requires entry of the ambient temperature and 
desired tempera ture of the strips. Then the strip number is 
entered and the program prompts for initial voltage drop and 
initial current of the strip. The program computes strip 
temperature using the formula: 




where C and S are taken from the calibration curves. 

Next the program calculates the voltage drop that would be 
expected at the desired strip temperature using: 



E 



D 



R^iWV 

J Ei R 



Cr -t ) + 
U D F' 



Ii 



D 



1/2 
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Where = required resistance, from calibration 

= Fluid temperature 
S = Slope of calibration curve 

The variac on the test rig is then adjusted to produce 
this voltage drop across the strip. The strip must be allowed 
to come to thermal equilibrium after an adjustment is made. 

The calibration curves assume thermal equilibrium and are mean- 
ingless when temperature transients are taking place. 

The program listed requires 204 program storage locations 
and 76 memory storage registers. The TI-59 must be reparti- 
tioned to obtain adequate storage registers by entering 8 (2nd) 

OP 17. This provides 319 program storage locations and 79 
data registers. All of the data listed must be entered to take 
advantage of the full capability of the program. 

After the program and stored data are loaded the ambient 
temperature is specified by entering the temperature and pres- 
sing key "B". The desired strip temperature can then be 
entered and the key "C" depressed. The program is now ready 
for use . 

In use the number of the desired strip is entered and key 
"A" depressed. The printer will then prompt "E?" and the initial 
voltage drop should be entered and the "R/S" key depressed. The 
program will echo the current on the printer and will proceed 
to print the current temperature and desired voltage drop. The 
program will then pause and is ready for the next strip number. 
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FIGURE 1 - SKETCH OF THE MODEL 
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FIGURE 2 - CROSS-SECTIONAL SKETCH OF THE MODEL 
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FIGURE 3 - PHOTOGRAPH OF THE MODEL 
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FIGURE 5 - PHOTOGRAPH OF THE WIND TUNNEL TEST 
SECTION AND INSTRUMENTATION 
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FIGURE 7 - PHOTOGRAPH OF- THE ROTATING SHUTTER 
VA1.VES 




FIGURE 8 - PHOTOGRAPH OF THE MODEL MOUNTED IN 
THE TUNNEL 
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CALIBRATION CIRCUIT 
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FIGURE 9 - SKETCH OF CALIBRATION AND WORKING 
CIRCUITS 
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FIGURE 10 - TYPICAL VELOCITY WAVEFORMS 
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FIGURE 11 - SKETCH OF TELEDELTOS PAPER MODEL 
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FIGURE 12 





- SKETCH OF THE LUMPED RESISTANCE 
MODEL 
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FREQUENCY = 5.00 CPS BLADES = 4.0 INCHES 

FREQUENCY NR = 0.464 x 10' 6 



RUN NK V - 9U DEGREE MODEL, HIGH FLOW RATE 
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VELCCITY = 139.4933B ET/S AMRLITUDL NK =0.115 

EREGUENCY = 22.00 CPS GLADES = 4.0 INCHES 

FREQUENCY NR = 1.166 x 10'£> 



RUN NR 10 - 90 DEGREE MODEL, LOW FLOW RATE 
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FREQUENCY = 50.00 CPS ULAUES = A.O INCHFS 
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FREQUENCY NR = 8.883 x 10" 
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FREQUENCY NR = 2.685 
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RUN NR 3 1 - 45 DEGREE MODEL, LOW FLOW RATE 
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KUN NR 32 - 45 DEGREE MODEL, HIGH FLOW RATE 
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FREQUENCY NR = 1.023 



RUN NR 33 - 45 DEGREE MODEL, LOW FLOW RATE 
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FREQUENCY NR = 1.757 



KUlM NK 34 - 45 DEGREE MODEL, LOW FLOW RAIL 
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RUN NR 35 - 45 DEGREE MODEL, LOW FLOW RATE 
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FREQUENCY NR = 3.452 x 10' 



RUG (MR 36 - 45 DEGREE MODEL, LOW FLOW RATE 
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RUN NR 37 - 45 DEGREE MODEL, LOW FLOW RATE 
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REYNOLDS NR = L07462.31250 STROUHAL NR = O.2H0 

VtLCCITY = 52.8C045 FT/S AMPLITUDE NR =0.650 
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ELLL1IY = 48.53128 FT/S AMPLITUDE NR =0.330 

KLCUEiMCY = 55.00 CPS bLADES = 6.0 INCHES 
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